Voltage-tunable singlet-triplet transition in lateral quantum dots 
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Results of calculations and high source-drain transport measurements are presented which demon- 
strate voltage-tunable entanglement of electron pairs in lateral quantum dots. At a fixed magnetic 
field, the application of a judiciously-chosen gate voltage alters the ground-state of an electron pair 
from an entagled spin singlet to a spin triplet. 
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I. INTRODUCTION 

Proposals for spin-base)i|-CHiantuni computation in 
a solid-state environmentld'ElBoO require efficient tech- 
niques for manipulating the entanglement of coupled 
qubits. In this paper, we demonstrate theoretically and 
verify experimentally that ground-state entaglement can 
be induced solely by applying a potential to the gates. 
This is possible because the gate voltage controls not 
only the chemical potential of the dot, but the shape of 
the confining potential as well. Consequently, the gate 
voltage can induce transitions in the dot containing a 
well-defined, constant number of particles. 

Early far-iiLfcared measurements on arrays of fewp. 
electron dots,l3iZl and transport through single devicesQ 
were focused on the tunability of the electron num- 
ber and, although transitions were observed, it was 
difficult to assign, for example, quantum numbers to 
these transitions. Later experiments using single-particle 
capacitanceB and magneto-tunnelingl3 spectroscopy fo- 
cused on the evolution of the ground states as a function 
of magnetic field and were able to distinguish features 
consistent with the two-electron singlet-triplet transition. 
High source-drain tunneling spectroscopy probes the ex- 
cited states as well as the ground state and therefore the 
singlet-triplet transition can be more clearly and unam- 
biguously observed. This has already been successfully 
applied to etchedj-wertical quantum dots with (spin) un- 
polarized leads OO In the lateral devices employed in the 
present study, the dot is formed within a 2-dimensional 
electron gas (2DEG), with the lateral confinement pro- 
duced electrostatically by voltages applied to gates lo- 
cated above_the 2DEG. Recent work employing a novel 
gate designtj has allowed the electron-number N to be 
tuned down to a single electron. 

There are at least two features unique in the lateral de- 
vices. First, since the leads are essentially 2DEG edges, 
applying a rather weak magnetic field — approximately 
0.4 T j« practice — is sufficient to produce spin-resolved 
edges.li3 Therefore, the tunneling rates into and out of the 
dots are significantly different for each species of spin. 
This spin-polarized injection (and detection) allows us 



to distinguish orbitaLjCfEccts. from spin effects in trans- 
port measurements.EjEdilj't3 Second, since the confine- 
ment potential is formed electrostatically by the various 
gate voltages, altering the shape — in particular the non- 
parabolicity — of the quantum dot can be accomplished 
while keeping the particle number fixed. It is impor- 
tant to note that the familiar singlet-triplet transition 
is not caused by the difference in Zeeman energy but 
rather by changes in the orbital part of the wavefunc- 
tion. In multiple-dot systems, in particular wiili regards 
to quantum-dot-based quantum computation,EI this ver- 
satility is of crucial importance. Quantum-state engi- 
neering of this sort is clearly observable in the experi- 
mental results we present below. 

From the theoreticaLaidn^umerical analyses of the in- 
teracting two-elecjtijpiJl3ll3llj problem, as well as higher 
electron numbersji3n3 demonstrated singlet-triplet tran- 
sitions in parabolic potentials using either a fixed or an 
A^-dependent harmonic frequency. These works focused 
on magnetic-field-induced transitions. In our experiment, 
the confining potential is a function of the continuous 
variables Vg (the gate voltages) and the confinement in 
our particular quantum dots deviate from parabolicity. 
These features, which are addressed in our theory, al- 
low a spin phase diagram to be constructed in the gate- 
voltage/magnetic- field plane and clearly indicate how the 
N = 2 singlet-triplet transition can be externally engi- 
neered at fixed magnetic field. In double-dot systems 
containing one electron apiece, similar transitions, and 
for similar reasmis, have been theoretically demonstrated 
in both lateralEilo and verticalEj devices. Our theory is 
applicable to these systems with only a few modifica- 
tions related to the orbital degrees of freedom — the spin 
physics are essentially equivalent. 

In the Coulomb-blockade regime, transport experi- 
ments probe the two-electron system either by adding an 
electron to the one-electron droplet, or by removing an 
electron from the three-electron droplet. Each case cor- 
responds to a distinct gate voltage, and in each case the 
ground and excited states can be probed by high source- 
drain spectroscopy, which directly reveals the singlet- 
triplet transition. Our theory and experiment show that, 
for these two different gate voltages, the transition be- 
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FIG. 1; SEM image of the gate geometry forming the quan- 
tum dot. This geometry enables a precisely known number 
of electrons (A*' = 0,1,2,... ,50) to be trappeda and pro- 
duces a quasi-parabolic confinement potential. Sweeping the 
plunger-gate voltage tunes both the shape and the chemical 
potential of the quantum dot. 



tween the entangled spin singlet |t|) — |J,t) and the spin 
triplet I J, J,) occurs at two different magnetic fields. 

The paper is organized as follows: Section II describes 
the experimental results, including the demonstration 
of a gate-voltage-induced singlet-triplet transition. Sec- 
contains a theoretical analysis. 



tion III contains a theoretical analysis, culminating in 
the spin phase diagram of the N = 2 interacting system 
in the gate- voltage/magnetic-field plane. Finally, Sec. IV 
contains a concluding discussion. 



II. EXPERIMENTAL RESULTS 

An SEM image of a device similar to the one used in 
our experiments is shown in Fig. |l|. This geometry allows 
us to controUably tune the number of trapped electrons 
in the 2DEG — 90 nm belom-the surface — from about fifty 
down to a single electron.Ej High source-drain transport 
measurements in the Coulomb-blockade regime were car- 
ried out in order to detect both the ground and excited 
states of the two-electron system. Standard low-power 
ac measurement techniques were used with a 10 /iV ex- 
citation voltage applied across the sample at a frequency 
of 23 Hz. An additional dc voltage was applied in or- 
der to obtain a high source-drain bias. The differential 
conductance dl/dVsd is measured directly in such a con- 
figuration and the relevant data is shown in Fig. |^ for a 
source-drain voltage of 350 fiY. In this figure we show the 
inverted grayscale for the N = 2 subspace as a function of 
magnetic field B and plunger gate voltage Vg. The neg- 
ative differential conductance, which is also tunable,E3 is 
related to the spin-polarized injection of electrons. 

For the lowest set of curves in Fig. ||, transport pro- 
ceeds through the addition and subtraction of a second 
electron from a one-electron droplet. At the lowest curve, 
transport is predominately through the ground state of 
the two-electron droplet (a spin singlet at low fields); be- 
ginning at the curve immediately above this one, trans- 
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FIG. 2: High source-drain transport spectroscopy of the two- 
electron droplet. The lower set of curves corresponds to fluc- 
tuations between A'^ = 1 and N = 2, while the upper set 
corresponds to fluctuations between N = 2 and A'^ = 3. Both 
sots of curves probe the same states of the two-electron droplet 
(ground plus first excited state) but at different gate voltages. 
The singlet-triplet transition is seen to occur at two different 
critical fields. The solid line marks the singlet-triplet ground- 
state boundary. This boundary can be traversed along the 
dashed line at B = 1 T solely by adjusting the (plunger) gate 
voltage. 



port through the first excited state (a spin triplet at low 
fields) is also allowed. Hence, the exchange constant can 
be directly obtained experimentally from these curves by 
suitably calibrating the parameters relating gate voltage 
to energy. The singlet-.tr.ipjet transition is clearly seen 
(c/. etched vertical dotsBB) at a field of B^i ~ 0.92 T. 
The upper set of curves corresponds to adding and re- 
moving a third electron from the two-electron system. 
After the third electron has left the dot, the resulting 
two-electron droplet can either be in the ground state or 
an excited state. (Transport through the ground state is 
the topmost curve.) Therefore, it should be possible to 
extract the same exchange constant from these curves as 
described above for the lower set of curves. Indeed, the 
singlet-triplet transition is again clearly seen, but now oc- 
curs at a field of Bc2 ~ 1-1 T. The singlet-triplet gaps for 
the two different cases are shown in Fig. ^ — the central 
experimental result of this work. 

Two important conclusions can be drawn from Fig. |^, 
each of which is verified in the subsequent sections. First, 
because the gaps do not close linearly, the confinement 
potential cannot be parabolic. Second, because the two 
curves do not fall on each other, the actual shape of the 
dot must be different for the two curves. This change can 
only be due to the gate voltage, and it therefore follows 
that the gates themselves can be used to tune through 
the singlet-triplet transition and hence tune the ground- 
state entanglement of the system. Figure ^ shows how 
this may be accomplished. The solid line demarcates the 
boundary between the singlet and triplet ground-state 
phases. At a fixed field of 1 T (marked on the figure as a 
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The first two terms in this equation are the oscilla- 
tor energies, with Q± — (a/w^ + (2wo)^ ± Wc)/2, lOc = 
eB /{•m*c) denoting the cyclotron frequency, and loq the 
parabolic confinement frequency. The final term is the 
Zeeman energy, where g^sB ~ 0.012hujc in GaAs. 

Neglecting environmental influences, the Hamiltonian 
of an isolated-quantum dot can be written in the Fock- 
Darwin basisEj as 



E Kf44'Cfe.'Q., (2) 



i,j,k,l, 

(7,(7 



FIG. 3: Singlet-triplet gap J — AE as a function of magnetic 
field for two different gate voltages, as described in the text. 
AEi and AE2 are also denoted in Fig. H. 



dashed line) the ground-state entanglement can be tuned 
to be either a singlet or a triplet solely by adjusting the 
gate voltage appropriately. 

In the following, we present the theoretical justification 
of the above statements. 



III. THEORETICAL RESULTS 

We begin this section with a description of the model 
we shall use throughout the paper. We shall work pri- 
marily in the 2D harmonic oscillator (Fock-Darwin) ba- 
sis, characterized by the two oscillator quantum numbers 
771, rt = 0,1,2,... and the spin quantum number a — 
±1/2. This is the diagonal basis of 2D electrons (taken 
to lie in the x-y plane) with charge — e and effective mass 
m* , moving in a uniform magnetic field B = (0,0, -B) 
oriented perpendicular to the 2DEG plane, and with a 
parabolic confinement potential V^ar = m*u!Q{x'^ +y'^)/2. 
The single-particle energy levels are gi'sf.en. hv the familiar 
Zeeman-split Fock-Darwin spectrumicJcj'EZl 



1 

^+2 



g^isBa. (1) 



where the Latin indices i, j,k,l are a composite denoting 
the two oscillator quantum numbers m and n. The op- 



erator c. 



creates a particle in the state 



with the z-component of angular momentum (m — n), 
and Cia- is the conjugate annihilation operator. Unless 
otherwise noted, all energies are expressed in units of 



the effective Rydberg Ry 



*e7(2r) (« 5.9 meV in 



GaAs), and all lengths in units of the effective Bohr ra- 
dius ao = l{m*e'^) (« 9.8 nm in GaAs). 

The diagonal one-body term — the first term in 
Eq. (^ — is just the Fock-Darwin energy coming from 
parabolic confinement; £^0- £mna is given in Eq. (|^). 
We consider the total confinement to be composed of a 
parabolic piece plus a non-parabolic piece; the parabolic 
piece, along with the kinetic energy, is incorporated into 
the diagonal term; the non-parabolic piece is represented 
by the off-diagonal one-body term — the second term in 
Eq. (^) — whose overall strength is governed by the di- 
mensionless parameter 7. This term directly affects the 
single-particle spectrum, and also significantly alters the 
singlet-triplet transition in the tw o-elec tron droplet. We 
discuss this term in detail in Sec. IIIB. 

Finally, the two-body term in Eq. (H) repre- 
sents interactions, where the matrix element V^^' — > 

(TOi,ni;m2,n2|e^/(e|fi - f2|)|'7i3, 'T'S; "^■4, "4) is the full 
Coulomb interaction in the 2D harmonic oscillator basis 
(e ia,the dielectric constant); an exact expression is given 
by@ 
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fc3=0 



V ^3 / V 



E 



fc4=0 



k2=0 



F k 



(3) 
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The energy scale of the Coulomb interaction is set by Eq = ^7r/2 e^/{e£o) where the hybrid length 
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i'^ = hc/{eBy^l + iujl/uj'^). Eo is the sum of all 
exchange energies in the lowest Landau level, i.e., 
'E'^=o{m,0;0,0\ey{e\ri-f2\)\m,0;0,0)=Eo, and, ad- 
ditionally, (0,0;0,0|eV(e|fi - f2|)|0,0;0,0) = Eq/V^. 
The Coulomb interaction conserves angular momentum 
R = ~ TT-i). This is enforced in Eq. (^) by the 

Kronecker delta function: Rl — (nii — '^■1) + — 7^2), 
Rr — {nT-a — "•3) + {^4 ~ "■4)- Also in Eq. (0), /c = 
(toi +m2 +^3 -l-n4) — (fci + ^2 + fca + ^4), and r(fc + 1/2) 
is the Gamma function. 

The dimensionless parameter a in Eq. (|^) controls the 
strength of the Coulomb interaction, with a = 1 rep- 
resenting "bare" Coulomb interactions. At long length 
scales, screening effects from the nearby metallic gates 
and leads decrease the strength of the Coulomb inter- 
action. At short length scales, the finite width of the 
2DEG layer also decreases the strength of interactions. 
Since Coulomb interactions are not the primary focus at 
present, we shall use the parameter < a < 1 to de- 
scribe the strength of Coulomb interactions, rather than 
considering a more sophisticated functional form. 

In order to determine which of the main results are 
specific to the details of the confinement and which are 
more general, we shall first consider the usual parabolic 
confinement, and subsequently investigate the particular 
deviations from parabolicity present in our device. 



A. Parabolic confinement 

In this section, we consider the case of pure parabolic 
confinement (7 = 0). We shall see that our central result 
of voltage-tuned entanglement is already present in this 
simple case. 



1. Non-interacting electrons 

The single-particle problem with parabolic confine- 
ment (7 = a = 0) yields the Fock-Darwin spectrum in 
Eq. ffl). This approximation should be most valid for 
the lowest-energy level of the one-electron droplet; in ad- 
dition to having no intra-dot Coulomb interactions, the 
zero-point energy should be smallest for the one-electron 
droplet. A comparison of the Fock-Darwin spectrum and 
experiment is shown in Fig. ^. The experimental points 
are the position of current peak as a function of magnetic 
field for the one-electron droplet. The parabolic confine- 
ment frequency loq was used as a fitting parameter, with 
ujQ — 1 meV being displayed in the figure. Although 
the data is well-fit for this value, we shall find a rather 
different situation for the two-electron droplet. 



2. Interacting electrons 

When Coulomb interactions are switched on (but with 
7 = for the moment), m and n are no longer good 
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FIG. 4: Comparison of experiment and theory for the one- 
electron droplet. The data was fit to Eq. (|l|), with ljo as a 
fitting parameter. 



quantum numbers, but, since circular symmetry is still 
manifest, the total angular momentum R is indeed con- 
served, as are total spin S and total Sz ■ The Hamiltonian 
can therefore be diagonalized in each (i?, 5, Sz) subspace 
separately. 

We have numerically diagonalized Eq. ^ with 7 = 
according to the following procedui^: We work in a 
fixed (i?, Sz) subspace — alternativelyJHj one may work in 
a fixed (i?, S, Sz) subspace, a particularly useful approach 
for larger particle numbers — and we use the 2D harmonic 
oscillator basis, with the Coulomb matrix elements given 
by Eq. (||). We truncate the infinite-dimensional Hilbert 
space by introducing a high-energy cutoff i^cutoff; For 
each Ai"-particle basis vector \miniai, . . . ,mNnN(yN) , we 
calculate its (a = 0) eigenenergy and discard it if this 
energy is greater than i?cutoff- We then numerically 
diagonalize the resulting finite-dimensional Hamilto- 
nian (mArriArcr at, .. . ,miniai\H\m'in[a[, . . . ,m'^n'j^a'^) 
(with finite a) to obtain both the eigenstates and the 
spectrum. We then keep repeating this process with a 
progressively larger i?cutoff until the eigenvalues converge 
to a constant value. 

For small magnetic fields (i? < 5 T, for ujq = 1 meV) 
and two electrons, convergence is reached rather quickly. 
For example, for a 112-dimensional Hilbert space, conver- 
gence to within 4% has been achieved for the lowest 65 
eigenstates for (TV, B, R, Sz) = (2, 0, 1, 0) (this spin value 
includes both singlet and triplet states), and to within 
0.5% for the lowest 49 eigenstates. The eigenvalues for 
R = -10 to 10, 5 = 0, and Sz = Q are shown in Fig. |. 

Experimentally, we have seen that the magneto- 
transport data of the one-electron droplet is very well de- 
scribed by parabolic confinement with a confinement fre- 
quency of = 1 meV. In the two-electron droplet, where 
Coulomb interactions are now relevant, the ground-state 
singlet-triplet transition is experimentally seen to occur 
at approximately Sc = 1 T. If we assume the confinement 
frequency remains constant at 1 meV, then Be = 1 T oc- 
curs for a ~ 0.2, and thus Coulomb interactions are sig- 
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FIG. 5: Eigenvalues of the interacting two-electron droplet 
with parabolic confinement. The plot is for B = T, ujo — 
1 meV, and for singlets and triplets with Sz = 0. 
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FIG. 6: Inset: Lowest-energy singlet |S, Sz) = |0, 0) (lower 
curve) and triplet 1) (upper curve) for the two-electron 
droplet with a = 1 and ojq = 1 meV. The singlet-triplet gap 
J is indicated. Main plot: Singlet-triplet gap as a function of 
magnetic field for the two-electron droplet with a = 1. The 
upper curve has wo = 1-5 meV and the lower has wo = 1 meV. 



nificantly reduced from their bare value. Alternatively, 
since the critical field scales with the ratio uJo/{aEo), ujq 
may increase, rather that a decrease, to give the same 
effect. This is shown in Fig. |6| where all curves are for 
the case a = 1. The main plot shows the evolution of 
the singlet-triplet gap J with magnetic field. The lower 
curve has t^o = 1 meV while the upper has = 1-5 meV. 
The inset shows the actual singlet-triplet crossing for 
ujQ = 1 meV. This simple model of parabolic confine- 
ment with full Coulomb interactions is clearly insufficient 
to quantitatively reproduce the experimental findings of 
Fig. ||. The linear closing of the gap in the theory — there 
are actually very slight deviations from linearity not dis- 
cernible in Fig. ^ — appears to be a feature of Coulomb 
interactions combined with parabolic confinement. The 
experimental curves of Fig. ^ are thus an indication that 
the confinement is non-parabolic. 



What this simple theory does capture, is that the 
critical-field i?c(wo) is controlled by the gate voltage via 
its influence on the conflning frequency uJo{Vg) for a fixed 
particle number N. Thus, already at this level, we see 
how the gate voltage, at fixed uniform magnetic field, can 
be used to tune through the singlet-triplet transition. 

The parabolic-confinement model is insufficient to re- 
produce both the critical field Be and the zero-field gap 
Jo = J{B — 0). Finite- width effects {a < 1) or an in- 
creasing luq serve to increase both Be and Jq. In the 
following section, we investigate the infiuence of non- 
parabolic confinement on the singlet-triplet gap. 



B. Non-parabolic confinement 

The confinement potential produced by the gate geom- 
etry shown in Fig. |l| exhibits only approximate circular 
symmetry and this explicit symmetry breaking can be 
clearly seen in experiment already at the two-electron 
level, as shown in Fig. ^. This deviation from parabolic- 
ity is included as the second term of our model in Eq. (||), 
and it is the infiuence of this term we investigate in this 
section. 

The first problem is to obtain a functional form for the 
confinement potential. Since the gate voltages are the 
primary contribution to the confinement potential, the 
simplest approach is to consider the electrostatic poten- 
tial in the 2DEG induced by the gate voltages. Defining 
the x-y plane to be the plane of the gates so that the po- 
tential V{r,z = 0) is experimentally pgiven (r = {x,y)), 
an analytic expression can be derivedE2l for the potential 
V{r, z) at an arbitrary point: 



V{r,z) 



dr' 
2^ 



y(r',0) 



r r' 



3/2 ■ 



(4) 



This equation yields the correct z ^ limit, as well as 
dV/dz for \z\ oo. The integration, performed at 
each point r in the 2DEG plane (z = 90 nm), yields the 
potential which laterally confines electrons. A contour 
plot of this confinement is shown in Fig. |^ We have 
neglected the contribution of any i3-dependant effects of 
the edge-states (i.e., the leads). In the Coulomb-blockade 
regime we are interested in here, the S-dependence of the 
lead states will primarily influence the tunneling rates 
into and out of the dot. That is to say, the amplitude of 
the current will be affected, rather than the spectrum of 
the dot. 

The confinement potential in Fig. ^ can be viewed as a 
sum of a parabolic dot and a parabolically-confined semi- 
circular wire of diameter D which intersects the quan- 
tum point contacts (seen as saddle-points in Fig. ^ and 
the center of the dot. These considerations lead to an 
analytic expression which very closely approximates the 
(numerically-derived) potential in Fig. |7. This potential 



is given by i/par+7-f^nopar, where i/par 



{ll2)m*ujl{x^ 
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FIG. 7; Contour plot (thin lines) of the two-dimensional con- 
finement potential produced by the gates (thick lines) located 
90 nm above the 2DEG. 



y^) is the usual parabolic confinement, and 



1 



TT ^ * 2 ( y 



(5) 



SV-SR {SVq is Hermitian), and where 

5Vo^3/3^{a''^d^ + 2aW^b + d^b^) (6b) 
+ (1 + 12f3^){aU + db) + (1 + 6/3^), 
2/3{aU^ + 2d + d^b), (6c) 
4(3^ {aU^ + d^b) + (1 - \2(3^)d^, (6d) 

(6e) 



5Vi 
5^3 = 2/3d^ 



and 



5Vi = P'd 



2j4 



(6f) 



with d = (a + fet), and f3 = D-\u^ + 4^^)-^/". 

Before going on to the main task of investigating the 
singlet-triplet transition in this non-parabolic potential, 
we first investigate how the Fock-Darwin spectrum is af- 
fected by these symmetry-breaking terms. 



is the non-parabolic piece of the confinement potential. 
In Eq. (^), the parabolic piece is incorporated into the 
diagonal one-body term, while the off-diagonal one-body 
term is the second-quantized version of Eq. (|^) with 



\m n 



^noparl 

The computational consequences of the explicit 
symmetry-breaking terms (7 7^ 0) are that the Hilbert- 
space truncation scheme must be altered to incorporate 
the mixing of different angular-momentum subspaces. 
One possibility for the interacting problem is to begin 
with the single-particle states (no longer simple Fock- 
Darwin states) and then solve the interacting problem in 
the exact single-particle basis. Another approach, indeed 
the one we present below, is to treat the parabolically- 
confined interacting problem (7 = 0, a 7^ 0) exactly, 
as was done in Sec. |lIA|, and, in this basis, treat the 



one-body symmetry-breaking terms. 

That Eq. (|) mixes different R subspaces can be ex- 
plicitly seen by re- writing the position oppsators in terms 
of the usual oscillator ladder operatorsO x = £o{a'^ + 
a + b^ + b)/V2 and y = -i4(at - a - &t + b)/V2, 



l)a 



t 



and b^ — 



where 

^)clyi_^_i ^Cm.n.a- Using these relations, we rewrite the 
second term of Eq. (0) as 



7 ^^Aa^^J'T 



7^0/4 



4 



^0 SR=0 



(5VsR + 5V^,j, 



(6a) 



where SVsr changes the angular momentum (m — n) of 
the single-particle state |mn) by an amount 5R, 5V^j^ = 



1. Non-interacting electrons 

We have been unable to find an exact analytic solution 
to the non-parabolically confined (7 ^ 0) single-particle 
{a — 0) problem, and we therefore employ a numeri- 
cal treatment. It is simplest to again work in the 2D 
harmonic oscillator (Fock-Darwin) basis. In the single- 
particle problem, Zeeman effects are rather trivial and 
shall therefore be neglected in the present discussion. 
We include a fixed number of Fock-Darwin states in our 
Hilbert space, diagonalize the Hamiltonian, and repeat 
with a larger number of Fock-Darwin states, progres- 
sively increasing the Hilbert-space dimension until con- 
vergence of the spectrum is attained for the lowest few 
levels. 

An example of the resulting spectrum is shown in 
Fig. |[ with the equivalent Fock-Darwin spectrum as an 
inset. The most dramatic effect of non-parabolicity is 
at low magnetic fields, where the shell structure is heav- 
ily renormalized, although its remnants are still observ- 
able. The plot was computed using the 1891 lowest- 
energy Fock-Darwin levels — corresponding, at zero field, 
to the first 61 shells in the parabolic case — and with 
ujQ — 1 meV, 7=1, and D = lOao. At much larger 7, 
the shell-splitting becomes so large that different shells 
overlap, and a Fock-Darwin description at low fields be- 
comes dubious. Apparent anti-crossings are also seen in 
Fig. 1^, whereas the Fock-Darwin spectrum contains only 
crossings. Another important feature is that the 1^ — 2 
line at moderate fields is still clearly visible, even for 
larger 7; beyond this point, field effects begin to play a 
more prominent role than non-parabolicity effects, whose 
presence in the spectra becomes concealed. Thus, we do 
not expect non-paijabolicity effects to play an appreciable 
role beyond = 2.E£l 
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FIG. 8: Main plot: Single-particle spectrum as a function 
of magnetic field for a non-parabolically confined dot with 
Wo = 1 meV, 7 = 1, and D = lOao. Remnants of shell- 
structure at zero field can be discerned. Inset: The same 
spectrum, but with 7 = 0, i.e., the Fock-Darwin spectrum. 



2. Interacting electrons 



In section III A, we treated the parabolicaUy-confined 



interacting case by first solving for the non-interacting 
(Fock-Darwin) case; these states were then used as the 
basis states in which the interacting problem was solved. 
Continuing the progression, we now use the exact inter- 
acting many-body states computed in Sec. Ill A as the 
basis states in which we treat the non-parabolic piece of 
the confinement potential. 

The method of truncating the Hilbert-space must 
again be chosen. As an example, Fig. || — for all angu- 
lar momenta since they are no longer conserved — may 
be considered the Hilbert-space for the particular case of 
LUo = 1 nieV, B = T, and S^, = 0. Two methods of 
truncating this Hilbert-space can be considered. First, 
the k lowest-energy states within each angular momen- 
tum channel may be chosen as the reduced Hilbert-space, 
with all higher-energy states discarded. For example, if 
k = 10 and there are 21 angular momentum channels (as 
in Fig. H), then the Hilbert space has 10 x 21 = 210 di- 
mensions. In this scheme, the cutoff energy is variable, 
but the number of states within each angular momentum 
channel is fixed. The second method is to employ a fixed 
energy cutoff i?cutoff and to allow a variable number of 
states within each angular momentum channel; all states 
above Scutoff, regardless of angular momentum, are dis- 
carded and all states below i?cutoff, regardless of angular 
momentum, are retained. In principle, it matters little 
which truncation scheme is used so long as each method, 
of course, converges to the same values. In practice, the 
second method, with a fixed cutoff energy, achieves con- 
vergence faster. 

Figure ^ shows results analogous to Fig. |^, but for the 
non-parabolic confinement discussed above. All curves 
shown in Fig. ^ have {a,j,D) = (1,3, 5ao). In the main 




0.4 0.6 0.8 1.0 1.2 
Magnetic Field [Tesla] 

FIG. 9: Main Plot: Singlet-triplet gap J as a function of mag- 
netic field for ujo ~ i- meV (lower curve) and luo ~ 1-5 meV 
(upper curve) for non-parabolically confined interacting elec- 
trons. Upper Inset: Lowest-energy singlet IS, ^z) = |0,0) 
(lower curve) and triplet |1, —1) (upper curve) as a function 
of magnetic field for cjq — 1 meV. Lower Inset: Shape of total 
confinement potential in real space. All curves in the main 
plot and in the insets are for a = 1, 7 = 3, and D — 5ao. This 
figure is the analog to Fig. ^ for non-parabolic confinement. 



plot, the gap J is plotted for luq — 1 meV (lower curve) 
and luq — 1.5 meV (upper curve). The upper inset shows 
the singlet-triplet transition for ujo = 1 meV, and the 
lower inset shows a sketch of the shape of the confinement 
potential for these parameter values. All these parame- 
ters, except for 7, are as in Fig. ^. The non-parabolic 
spectrum differs significantly from the parabolic spec- 
trum, particularly at small field, where the triplet has a 
much weaker field dependence in the present case relative 
to the parabolic-confinement case. This is the behavior 
also seen in experiment. (See Fig. ^j.) 

In general, the non-parabolic model yields results much 
closer to experiment than parabolic confinement, and it is 
clear that the confinement potential in experiment is not 
parabolic. The available phase space — with variation in 
wq, a, 7, and D — is rather large and so an optimal fit has 
not been performed. Nevertheless, the theoretical curves 
in Fig. ^ are in qualitative agreement with the experimen- 
tal curves in Fig. ^. Future analyses of high source-drain 
spectroscopy with higher electron numbers will produce 
additional constraints which will more meaningfully re- 
duce the range of these parameters. 

The central point is that the particular critical field 
Be obtained is itself a function of the confining-potential 
parameters, loq, 7, and D. Since each of these parameters 
is controlled by the voltage on the gates shown in Fig. [|, 
it follows that the gates themselves can be used to tune 
through the singlet-triplet transition and hence tune the 
ground-state ent ang lement of the system. This is shown 
explicitly in Fig. |lO| where we plot the spin phase diagram 
in the luq-Bc plane. In this figure, the singlet-triplet gap 
J is computed for various values of luq (tunable by the 
gate voltage); the critical field Be is then extracted from 
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FIG. 10: Spin phase diagram for the two-electron system 
in the loq-Bc plane. This plot has a = 1, 7 = 3, and 
D = 5ao, and loq is controlled by the gate voltages. The 
usual field-induced singlet-triplet transition and the voltage- 
induced transitions are indicated. 

the solution to J{u!q,Bc) = 0. In practice, the range of 
luq which is experimentally accessible is delimited by the 
TV = 3 and N = 1 subspaces as the gate voltage is swept. 
However, if the accessible range is sufficient to be seen by 
experiment, a voltage-tuned singlet-triplet transition at 
fixed B is achievable. As shown in Sec. II, this is indeed 
the case for the present experiment. 

IV. DISCUSSION 

The fact that the gate voltage (or set of gate voltages 
in the present case) controls not only the chemical po- 
tential of the dot but also the shape of the dot has been 
exploited in the present work to tunc the ground-state 
entanglement of an electron pair in a lateral quantum 
dot. The experimental evidence clearly shows that the 
two-electron singlet-triplet transition occurs at a criti- 
cal field which depends on the gate voltage. The con- 



finement potential is not parabolic (nor elliptical) and 
this allows great flexibility in changing the shape of the 
dot while simultaneously keeping the number of con- 
fined electrons fixed. However, the experimental demon- 
stration of voltage-tuned entanglement is not dependent 
upon the precise shape of the potential, and should be 
achievable in a wide range of potential shapes. The main 
requirement is that the gate voltages appreciably change 
the shape of the potential while the particle number re- 
mains constant. The present is a modest step towards the 
construction of a quantum gate in a solid-state environ- 
ment; it does not, for example, demonstrate the unitary 
evolution of the system between singlet and triplet states. 

Because the tunneling barriers into and out of the dot 
are large, only the leading order contribution to the tun- 
neling current can be seen in our experiment. In prin- 
ciple, however, the tunnel barriers can be reduced while 
still remaining in the Coulomb blockade regimejn order 
to measure an appreciable cotunneling current .til In this 
way, the singlet-triplet transition may be experimentally 
probed throughout the N — 2 subspace, and thus give a 
more stringent test of theory. 

In a double-dot system with one trapped electron 
apiece, essentially the same spin physics occurs, and thus 
a voltage-tuned ground-state transition should also oc- 
cur along with the consequent implications for quantum 
computing. Alternatively, the current work may be spec- 
ulatively viewed as a possible gate-controlled singlc-qubit 
operation, where the single coded qubit exists in a single 
quantum dot containing two (or more) electrons. This 
conjecture will be more fully developed in a future pub- 
lication. 
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